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^ ■ Abstract 

.W. 

We investigate the phase structure of three-dimensional quantum gravity coupled to 
an Ising spin system by means of numerical simulations. The quantum gravity part 
is modelled by the summation over random simplicial manifolds, and the Ising spins 
are located in the center of the tetrahedra, which constitute the building blocks of 
the piecewise linear manifold. We find that the coupling between spin and geometry 
is weak away from the critical point of the Ising model. At the critical point there 
is clear coupling, which however does not seem to change the first order transition 
between the "hot" and "cold" phase of three dimensional simplicial quantum gravity 
observed earlier. 
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1 Introduction 



In the recent years there has been a remarkable progress in our understanding of 
two-dimensional quantum gravity coupled to matter. Although it is still somewhat 
unclear what kind of theory we deal with when we consider pure 2d quantum gravity 
(how is the "average" intrinsic geometry to be characterized etc.), we get definite 
answers, if we ask for the modification of critical exponents of conformal field theories 
when they are coupled covariantly to 2d quantum gravity. From the point of view of 
string theory it is important to study the two-dimensional models, but if we want to 
consider possible theories of quantum gravity outside the context of string theories 
(and there are good reasons for taking such a point of view) they are to be considered 
only as toy models for higher dimensional gravity. In this letter we attempt to take 
a first step in the direction of the study of higher dimensional gravity coupled to 
matter. We will confine ourselves to three dimensions. There are several reasons 
for this. Firstly we do not have available the same powerful analytical methods as 
in two dimensions and we have to rely on numerical simulations. These are quite 
a lot easier in three dimensions than in four dimensions. In addition we expect 
three dimensional gravity to be placed in between the solvable two-dimensional case 
and four-dimensional gravity in the following sense: Three-dimensional (classical) 
gravity has no dynamical degrees of freedom by itself, as is also the case for two- 
dimensional gravity. But as for two-dimensional gravity this does not imply that it 
cannot couple in a non-trivial way to the matter fields. On the other hand the theory 
shares with the four- dimensional theory the feature that the Einstein-Hilbert term 
in the action is non-trivial, not renormalizable and (if not regularized) unbounded 
from below. From this point of view we might get important hints which can be used 
in four dimensions from a theory, where we can hopefully use some of the intuition 
we have now gathered in the two-dimensional studies. Furthermore it is much easier 
to find non-trivial statistical models in three- than in four dimensions. In fact the 
simplest example is the Ising model which we are going to use. 



Three-dimensional quantum gravity modelled on random triangulations was intro- 
duced in refs. Jj], [| |J and numerical simulations first performed in refs. @, ||, ||, [7], |j . 
Let us briefly describe the model and summarize the results. 

The continuum action of Einstein-Hilbert gravity in three dimensions can be 
written as 



2 The model 
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In the dynamical triangulated approach the functional integral over metrics is re- 
placed by a summation over all possible triangulations. Here, as in two dimensions, 
it is important that we restrict ourselves to a fixed topology (which we take to 
be that of S 3 ). The building blocks which we glue together in order to form the 
piecewise linear manifolds are regular tetrahedra. Consequently we will get (see for 
instance for a careful discussion) 

Jd 3 ^ ~ N 3 (2) 

J d 3 ^R ~ cm-QNs (3) 

where N Q , N 3 denote the number of vertices, links, triangles and tetrahedra which 
constitute the simplicial manifold. The number c = 2ir/ arccos(l/3) is chosen such 
that R = corresponds to flat three-dimensional space. This means that the dis- 
cretized action can be written as 

S[T] = k 3 N 3 - fcjJVi (4) 

and the recipe for going from the continuum functional to the discretized one will 
be: 

/PWe-^-E^. (5) 

The result of the numerical simulations is as follows: For small fci (i.e. for large bare 
gravitational coupling constant) the system seems to be in a phase of large Hausdorff 
dimension (denoted the "hot" phase). Even for large systems consisting of 28.000 
tetrahedra the linear extension (the average geodesic distance) of the system is small 
and increases only slowly with volume. For k\ ~ 4.0 there is a phase transition to 
a state where the system is quite extended. The extension grows linear with the 
volume, showing that the Hausdorff dimension is one. This phase is probably a 
lattice artifact which signals the dominance of the conformal mode of gravity for a 
small gravitational constant. The transition between the two phases is a first order 
transition since pronounced hysteresis is observed. We conclude that if we restrict 
ourselves to actions of the form (|4j) it seems not possible to define a continuum limit 
of the lattice model in the usual sense, with a divergent correlation length. This 
is not an undesirable situation since the presence of a divergent correlation length 
would force us to identify at least one massless field in three-dimensional quantum 
gravity. But we know such fields are not present in classical three-dimensional 
gravity. 

Let us couple the above defined model to Ising spins. The coupling is done in 
the same way as for discretized two-dimensional quantum gravity, where we know 
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the dynamical triangulated model coupled to Ising spins leads to critical exponents 
which agree with the ones calculated using continuum formalism. To each tetra- 
hedron i we associate an Ising spin cij and the partition function for the combined 
system is 

Z(p, h, k 3 ) = ]T e~ k ^ E e k ^e^<^> {& ^- 1] . (6) 

7V 3 T~N 3 [a] 

In this formula T ~ N 3 signifies the summation over all piecewise linear mani- 
folds which can be formed by gluing N 3 (regular) tetrahedra together such that the 
topology is that of S 3 . means the summation over all spin configurations while 
Yj<i,j> stands for the summation over all neighbour pairs of tetrahedra. One annoy- 
ing aspect of the above formalism is that we are forced to perform a grand canonical 
simulation where N 3 is not fixed. The reason is that we (contrary to two dimensions) 
have no ergodic updating algorithm which preserves the volume N 3 . In practise it 
is however possible to perform the measurements at a fixed N 3 and the important 
coupling constants will then be (3 and k\. We refer to |7j for a detailed discussion. 
The spin updating is performed by the single cluster variant of the Swendsen-Wang 
algorithm developed by Wolff M. The cluster updating algorithms have been suc- 



cessfully applied to the Ising model coupled to 2d gravity [10], [IT], 0, [L3| and to the 
ordinary three dimensional Ising model |l4j . 



3 Numerical results 

As mentioned above three-dimensional simplicial quantum gravity has two phases 
depending on the value of k\. The first statement we can make is that this is 
unchanged by the coupling to Ising spin. 

In the "hot" phase [k\ < 4.0) where the Hausdorff dimension is large the numer- 
ical value of the magnetization 

(7) 

is shown in fig.l as a function of (3. We see a clear signal indicating a phase transition 
from a disordered phase (small (3) where \a\ ~ to an ordered phase (large f3) where 
| a | w 1. The transition becomes sharper with increased volume N 3 and seems to be 
a second order transition. This situation is contrasted by the magnetization curve 
in the "cold" phase shown in fig. 2. Here is only a gradual cross over to |er| ~ 1 for 
large (3, and the cross over is weakened for increased volume N 3 . The situation is 
precisely as one would expect in the case of a one-dimensional system where there is 
no spontaneous magnetization. We conclude that the Hausdorff dimension d H « 1 



(7 
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measured in the pure gravity case seems to reflect correctly the dimension relevant 
for coupling to matter. 

In the k\ — [3 plane we have the phase-diagram shown in fig.3. If (3 is away from 
the critical value /3 c (ki) (which has only a weak dependence on fci) the coupling 
between the fluctuations in geometry and spin seems weak and of course it vanishes 
in the limits (3 — > oo and (3 — > 0. In these limits we therefore have a strong first order 
transition between the "hot" and the "cold" phase of three-dimensional quantum 
gravity, precisely as is the case in the absence of spins ||. In the "hot" phase, where 
the Ising system has a second order transition, we have seen an increased coupling 
between geometry and spins when we approach the critical (3 c (ki). This is shown in 
fig.4 where we plot the average curvature (R) as a function of f3. A clear peak is seen 
at j3 c . This enhanced coupling between geometry and spins at the critical point is 
qualitatively in agreement with the 2d results, where we have a change in the string 
susceptibility exponent ^string (not to be confused with the magnetic susceptibility 
exponent 7 ma9 ) from the pure gravity value — 1/2 to —1/3, precisely when (3 = (3 C . 
Unfortunately it is not clear that the entropy exponent analogous to "^string exists 
in the hot phase of three-dimensional quantum gravity (@, H [7|, H) so we have no 
obvious exponent with which we can compare the effect of the spin coupling, but the 
enhanced coupling between spin and geometry leaves open the possibility that the 
transition between the "hot" and "cold" phase changes from a first order to a second 
order transition. We have looked for hysteresis when changing k\ and adjusting (3 
to the critical value f3 c (ki). While the hysteresis is indeed weaker when measured 
this way, we still see a clear hysteresis (fig. 5) and we conclude that there is never a 
second order transition in geometry. 

Let us make the following remark concerning the determination of the phase 
diagram shown in fig.3: Due to the strong hysteresis it is somewhat ambiguous. We 
have used the following procedure: Well inside the "hot" phase the system follows a 
unique path when changing k\ and keeping (3 fixed as illustrated in fig. 5. The precise 
location depends on the value of (3. We have extrapolated these paths until they 
intersect the parts of the hysteresis curves which correspond to the "cold" phase. 

In fig. 6 we have shown the spin-spin correlation as a function of geodesic distance. 
To be precise there are two obvious candidates for geodesic distances (see |L5| for a 
discussion in the context of four-dimensional simplicial quantum gravity). We can 
define the geodesic distance d\ between two vertices as the length of the shortest 
path along links connecting the two vertices. Alternatively we could have defined 
the geodesic distance d 2 between two tetrahedra as the length of the shortest path 
connecting the two tetrahedra, moving from center to center in neighbour tetrahedra 



5 



which have a triangle in common. The first definition is clearly much closer to the 
"correct" definition obtained by considering the manifolds as piecewise linear, with 
the curvature attached to the links. The other definition corresponds to moving 
along links in the dual graph, which is a 4 graph. For a single manifold the two 
distances can differ a lot, but when an ensemble average is taken it seems as if 
one can consider them as proportional. A similar result is true in four- dimensional 
gravity (JOI). Clearly d 2 is most convenient for our purpose and starting from a 
given tetrahedron i we define the volume Vi{r) inside a ball of d 2 -geodesic radius r 
around %q as the number of tetrahedra within this distance. Further the differential 
volume is dVzfr) = V^r) — V^(r — 1). We can now define a spin-spin correlation 
function as 

An alternative correlation function would be 

G(r) = ( £ amX (9) 

which is related to the magnetic susceptibility x(P) by 

J2G(r)= X ((3)-\(3-(3 c (k 1 )r for /? - (10) 

r 

In fig. 6 we have shown g(r) for two different values of j3. In principle one can 
extract the mass gap m(/3) from the exponential fall off of g(r) and in this way 
determine the critical exponent v defined by m(/3) — \/3 — /3 C \ U . This seems however 
difficult to do in a reliable way, in accordance with the experience in two-dimensional 
quantum gravity coupled to Ising spins, and it is maybe understandable if one keeps 
in mind that not only is a precise determination of f3 c needed in order to extract 
v. In addition our data are folded into the distribution dVi{r) which determines 
the Hausdorff dimension, a quantity which by itself is very difficult to measure in a 
reliable way. 

In the same way we can construct x{0) from G(r), but it does not lead to a precise 
determination of the critical exponent 7 (again in agreement with the experience 
from two dimensions). 

4 Discussion 

We have shown that the phase structure of three-dimensional simplicial quantum 
gravity, as described in [Q, is not modified by the presence of matter, at least in the 
simplest case of coupling to Ising spins. In the so-called "hot" phase the Ising spin 
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system has a second order transition, while it has no transition in the "cold" phase 
in agreement with the effective one-dimensional nature of this phase. The existence 
of two phases in three-dimensional gravity is caused by the Einstein- Hilbert term in 
the action. In two dimensions this term is absent (for a fixed topology) and we have 
only one phase of pure gravity. This phase seems to have most in common with the 
"hot" phase of three-dimensional gravity and it is natural to expect that the critical 
properties of the matter theories covariantly coupled to 3d gravity could be changed 
in a non-trivial way when we are in this "hot" phase, simply by analogy with the 
two-dimensional models. 

One way to investigate the possible non-trivial scaling of the Ising model in 
the "hot" phase is by means of finite size scaling. In two-dimensional gravity this 
approach seems to work somewhat better than the direct attempts to measure v 
and 7 mentioned above. The disadvantage of the method is that it only gives us 
certain combinations of the exponents. 

If a given thermodynamic function F has a critical behaviour 

F(fl~(P-/3e)- (11) 

one expects in ordinary flat space a finite size dependence of the form 

F(/3,L)=L*f(\j3-/3 e \L 1 ' v ) (12) 

where L denotes the linear size of the system and the exponent v is determined by 
the divergence of the correlation length £(/3): 

m~\p-Pc\-'. (13) 

By a measurement of F(/3 C ,L) ~ /(0)L* as a function of L we can determine the 
combination x/u, while measurements away from (3 C would give us x directly for 
sufficiently large L. 

If we want to use these formulae for systems coupled to quantum gravity we 
must identify the divergent correlation length £ in terms of geodesic distances, as 
was already discussed in the last section. We further have to identify the linear 
extension L. If the system has a finite Hausdorff dimension du it is tempting to 
define 

L~Nl /dH . (14) 



These ideas have been used with some success in 2d-quantum gravity [11], [12], [13| , 
but could be spoiled if the Hausdorff dimension is infinite. Since the dn from (very 
tentative) direct measurements of V^(r) seems large in the "hot" phase it is tempting 
to conjecture that the critical exponents of the Ising model in the "hot" phase take 
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their mean-field values. A reliable determination of the various exponents along the 
lines discussed above requires a considerable amount of computer time, since it is 
already quite demanding in two dimensions, but we hope to be able to address the 
question in a future publication. 



Note added: While completing this article we received a paper by Baillie [IB], who 
investigates the same system. There seems to be little overlap with our work since 
he did not explore the phase structure in k\ — /3-plane. In fact this is not possible 
with the size of systems he uses. 
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Figure Captions 

Fig.l The magnetization \a\ (denned by (|7|)) as a function of ft in the "hot" phase 
for iV 3 = 4000 (triangles) and iV 3 = 10000 (circles). 

Fig. 2 The magnetization \a\ (defined by (0)) as a function of ft in the "cold" phase 
for iV 3 = 4000 (triangles) and iV 3 = 10000 (circles). 

Fig. 3 The phase diagram for 3d quantum gravity coupled to matter. Filled circles 
are results obtained for N 3 = 10000. 

Fig. 4 The average curvature (R) as a function of ft for iV 3 = 4000 (circles) and 
iV 3 = 10000 (squares). The position of the peak coincides with the value of ft c 
determined from the magnetization curve. 

Fig. 5 The hysteresis curve for pure gravity (triangles) and in the case where the 
Ising spin system is critical i.e. where it couples in a maximal way to gravity 
(circles). N 3 = 10000. 

Fig. 6 The spin-spin correlation function g(r) (defined by (|8|)) as a function of the 
geodesic distance r for ft = 0.5 (full drawn curve) and ft = 0.8 (dotted curve), 
k\ = 3.7 and N 3 = 4000. The best estimate of the critical value of ft is: ft c = 
0.85 
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